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We investigate a dissipative, deterministic ratchet model in the over- 
damped regime driven by a rectangular force. Extensive numerical calcula- 
tions are presented in a diagram depicting the drift velocity as a function of 
a wide range of the driving parameter values. We also present some theoret- 
ical considerations which explain some features of the mentioned diagram. 
In particular, we proof the existence of regions in the driving parameter 
space with bounded particle motion possessing zero current. Moreover, 
we present an explicit analytical expression for the drift velocity in the 
adiabatic limit. 

PACS numbers: 05.60-k, 05.45.Pq, 05.45.Ac, 05.45.Xt 

1. Introduction and model set-up 

Directed current in ratchet systems have received much attention over 
recent years JJ^OII^IS]. One of the reasons to study these type of systems 
is motivated by the attempt to understand the physical mechanism of mo- 
tion of molecular motors in biological systems and to investigate its role 
in the design of new material properties [SI 15 • The dynamics of a particle 
moving in a periodic potential under the action of an applied time periodic 
term is rather complex and rich, displaying typically even chaotic behavior 



* This work is dedicated to Professor Peter Talkner on the occasion of his 60th birthday. 



(1) 



2 



template'vS printed on February 6, 2008 



[3 IHl El ^1 • A feature of particular interest is the emergence of directed 
currents in the system response to an external time-periodic driving with 
zero time-averaged value. Though most of the works about ratchet systems 
consider the presence of noise, this phenomenon may also arise in deter- 
ministic systems, both in the overdamped ^2 ^1 and underdamped 
[71 IS] regimes. Moreover, the phenomenon of anticipated synchronization 
occurring in inertial deterministic ratchets have been studied recently |14j . 

In particular, in the overdamped regime, the one-dimensional particle 
dynamics x{t) usually considered is governed by a first order differential 
equation of the type 

x{t) = -U'[x{t)]+F{t), (1) 

where the dot and the prime denote time and spatial derivatives, respec- 
tively, U{x) is a periodic potential with spatial period A [i.e., C/(x + A) = 
U{x)], and F{t) is a time-periodic driving force with period T [i.e., F{t + 
T) = F{t)]. In this type of systems, the current is defined as the average 
velocity 

x(t)-x(0) 

V = lim — —. (2) 

As shown in Ref. jllj co-existing attractors can exist for large driving 
strengths, which however are not current-carrying. A finite current u, pos- 
sessing an unbounded x(i)-trajectory, is consequently independent of the 
initial condition x(0). It can also be shown that a directed current {v ^ 0) is 
only possible if at least one of the following symmetries is broken pilHISlIT^: 

3x0 £ ^ such that U{xq — x) = U{xq + x) Vx G SR, (3) 

F{t + r/2) = -F{t) Vt G 3f?. (4) 

Our main interest in this paper is to gain a deeper insight into these type 
of systems exploring some quantitative and qualitative aspects of its very 
rich dynamics. Specifically, we will consider the same sawtooth potential as 
in Refs. dH El 



sin(27rx) + — sin(47rx) 



(5) 



which has a spatial period A = 1. In addition, rather than using a sinusoidal 
driving force F{t)^ in this work we will consider a multifrequency time- 
periodic force given by 

T?(+\ / ^ for < t < I 
^(*)= -A for %<t<T, 
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Fig. 1. Plots of the instantaneous velocities —U'{x)+A (solid line) and —U'{x) — A 
(dashed line) as a function of x for a driving force with A — 0.5. (All quantities in 
dimensionless units). 

with A being a constant. Since the potential © breaks the spatial symmetry 
(jni), a directed current is possible, even though the time symmetry (jlj is 
fulfilled. 

The paper is organized as follows. In Section [2 we perform a detailed 
numerical study of Eq. for a wide range of the driving parameter values 
A and uj = 2tx jT . In particular, we present a colored phase-diagram v vs. 
A and w, which is inspired by the figure provided by Prof. Peter Talkner 
and collaborators in jl6j . In order to explain some features of this diagram 
we propose in Section |51 some simple theoretical considerations. Finally, in 
the last section, we summarize our findings. 

2. Dynamical regimes and numerical evaluation of the current 

Regardless of the value of the driving frequency, there exist two critical 
amplitudes = 3/4 and Al^ =3/2 separating three different regimes. For 
amplitudes A G [0, A*], both potential states, \J{x) + Ax and \J{x) — Ax, 
possess a periodic array of equilibrium points (see Fig. Since the inertial 
term mx{t) is absent, the particle cannot cross these equilibrium points, and 
consequently, it remains trapped between them, leading to a zero current. 
For A £ {Al,A2], the potential state U{x) + Ax retains its equilibrium 
points, while U{x) — Ax does not have any, that allowing non-bounding 
motion in the positive direction. Finally, when A > A2, neither of the 
potential states possesses equilibrium points, and the particle motion is not 
bounded in either direction. 

In order to go further in our analysis, we have resorted to a numerical 
treatment of our model using the freely available integrator RKSUITE . 
Wide regions of parameter space have been explored. In Fig.|21we have used 
a color code to represent the drift velocity v for a rectangular driving force 
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Fig. 2. The "Talkner" -current phase diagram. The drift velocity as a function of 
the driving parameters A and uj is represented using a color density plot. Black 
color has been used for velocities less than 10~^. (All quantities in dimensionless 
units) 

of varying amplitude A and frequency w. The regions in black correspond to 
bounded, time-periodic particle motion with zero current. More precisely, 
the velocity in those regions is smaller than 10~^. The diagram has a rich 
and interesting structure. Regions of zero drift velocity are intermingled 
with regions with finite values, giving rise to a finger-shaped structure. No- 
tice that for the parameter values considered, the drift velocity is always 
positive or zero, showing the largest magnitude in the intermediate region 
A £ {Al, A2]. This last feature could be understood by taking into account 
that in this regime the motion in the positive direction is never compensated 
by motion in the opposite direction. 

In Fig. |31 we present a section of the diagram for a driving amplitude in 
the intermediate regime A = lA. A series of peaks corresponding to the 
fingers in Fig. |21 are observed. Figure ^ shows a representative section for 
A > A2. By contrast with Fig. |51 it presents an intermediate gap of particle 
localization. If we further increase the value of A, more intermediate gaps 
appear, as can be seen in the phase diagram. 
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Fig. 3. Drift velocity v (solid line) as a function of the driving frequency uj for 
A = 1.4 < Ag. The dashed line shows the adiabatic value. The arrows depicts the 
zero- velocity bands given by the theory in the text. (All quantities in dimensionless 
units.) 




Fig. 4. The same as in Fig. El but for A = 2> 



3. Some theoretical results 



Even though the nonlinearity of the system precludes a complete and 
detailed analytical solution of the problem, it is possible to explain some 
features of the phase diagram by simple considerations. 
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Fig. 5. Determination of zero- velocity bands with x{t + T) = x{t). Solid line 
depicts G{x), defined in Eq. (|f 0|l . as a function of x for the same value of A as in 
Fig.EI 

3. 1 . Proof of existence of regions displaying particle localization 

Due to the truncation implicit in any numerical calculation, the simula- 
tions reported above are not able to distinguish between a situation of exact 
particle localization or a very small drift, lower than the tolerance we chose 
to determine the black regions in Fig|2 Furthermore, a simple explanation 
of the existence of such regions in the non-trivial regime A> A\ would be 
desirable. 

Let us assume that the particle starts at xq at the beginning of a driving 
period. It then moves under a force —U'{x) + A for half a period until it 
reaches the position xi > xq. Then the driving force switches sign so the 
total force on the particle is now —U'{x) — A, which makes it (in general) 
move backwards up to a position x'q after another T/2. If Xq = xq, the 
particle has returned to the initial position, and consequently, the process is 
repeated successively, leading to a drift velocity strictly equal to zero. For 
this situation to happen the following equations must hold 




T 




dx 



(7) 



2 



U'{x) - A' 



(8) 



Subtracting both equations we arrive to the condition 



G(xo) = 



(9) 
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where 

Therefore, we can determine the set of pairs {xq,xi) with zero current for 
a given driving strength A by plotting the function G{x) vs. x (see Fig[SJ. 
The intersection of a horizontal line with G{x) in this plot provides the 
possible values of the pair (xq, xi). The period T associated with the pair is 
then given by Eq. ((T)) or (jSJ. Since the drift velocity v is independent of the 
initial conditions [15) . it would be exactly zero for those driving parameter 
values A and T. 

In Fig. IHl we plot G{x) for the same value oi A = 2{> A2) as in Fig. |1| 
A horizontal line crosses G{x) at the points A, B, C, and D, providing the 
coordinates xa,xb,xc, and xd- If we choose xq = xa, and xi as any of the 
other points, we obtain three pairs of points with driving periods Tab,Tac, 
and Tad- Since the integrand in Eq. ((T)) is always positive, the further away 
xi is from xq, the larger the period, which implies Tab < Tag < Tad- 

Before we proceed in analyzing the plot in greater detail, let us study 
the symmetries of G{x). Since it is an integral of a space-periodic function 
[see Eq. (|10|) ] we have 

G{x + l) = G{x) + cp, (11) 

where (j) = In addition, taking into account that U'{x) is an even 

function, necessarily 

G{-x) = -G{x). (12) 

Property (|llj) leads to the fact that we only need to vary xq within a spatial 
interval of length unity, whereas property (|11() combined with H12() implies 
that G{x) has an inversion center at x = n, with n any integer. The driving 
period T, when viewed as a function of xi or xq [see Eq. (jTj)], also obeys 
these properties. 

Because of these symmetry properties, choosing the point A as xq in 
Fig.lSJis equivalent to choosing any of the points A', A", and A'", in the sense 
that they lead to the same driving periods. Let us consider a horizontal line 
between points A and A'. As this line gets closer to A, two of the intersection 
points approach each other, and consequently, the driving period associated 
to them tends to zero. Therefore, T = (w = 00) corresponds to particle 
localization. Moving up the line continuously we obtain a whole band of 
driving periods starting from zero up to a maximum value given by the 
pair {xbi-,xai) when the line crosses A'. This pair is equivalent to {xa-,xb)- 
Two more intermediate periods in the band can be obtained from the pairs 
{xbiXc) and {xc-,xd) in the figure. 

The next value of the driving period obtained from Fig. [3] is the one 
given by the pair {xa, xq)- This pair marks the start of a second band that 
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Fig. 6. Drift velocity v (solid line) as a function of the driving amplitude A in the 
adiabatic limit (tj — > 0). The curve is indistinguishable from the numerical results 
obtained with uj = 0.01. (All quantities in dimensionless units.) 



ends at the maximum value given by the pair {xa,xd)- Moving up the hne 
up to A" gives ah the intermediate values in the band. 

Clearly, the first band is associated with movement within a spatial 
interval less than unity (i.e., < xi — xq < 1), whereas the second band 
is related to displacements two times that distance (1 < xi — xq < 2). In 
Fig.|l]we have indicated with arrows the calculated frequency bands. It can 
be seen that they do not cover the entire region with numerically evaluated 
zero velocity. This is due to the fact that the above discussed mechanism 
is just the simplest one leading to zero current. The drift velocity can also 
vanish because the particle returns to its initial position after two or more 
driving periods, instead of after the first one. 

The same analysis can be carried out for a subthreshold driving A G 
(^^,^2]- III this case G{x) presents singularities at the equilibrium points, 
which prevents the particle from crossing those points when F{t) = —A. 
This leads to a single band which is near the zero-current region observed 
in the simulations, as shown in Fig. El 



3.2. Adiabatic limit 

In this section we will obtain an analytical expression for the current in 
the adiabatic limit a; — > 

v^{A) = hmv{A,uj). (13) 
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Specifically, as proved later on, Uad(^) is given by the average value 

vUA) = ]^[v^{A) + v.{A)], (14) 

where Vj^{A) and V-{A) are, respectively, the velocities in the presence of 
the static forces —U'[x) + A and —U'{x) — A. Obviously, Vj^{A) = for 
A G [0,A*] and v^{A) = for A G [0,^3], as the particle ends up being 
trapped by an equilibrium point. For A> A\, 

MA) = (15) 

where 



^ dx ^ 2^2A + ^-3 + AA{l + A) 

^ Jo A-U'{x) ^{-1 + 2 A) {3 + 2 A) {-3 + 4 A) ^ ' 

is the time taken for the particle to travel a distance equal to 1 (the spatial 
period) in the presence of the static force —U'{x) + A. Analogously, for 
A > A*2, 

MA) = (17) 

where 



r (A) - _ V4A + 2^3 + 41 + ^ AA - 2^3 + 41 

'~Jo A + U'{x)~ v/(-3 + 2A)(l + 2^)(3 + 4A) ^ ' 

is defined as t^{A) but replacing —U'{x) + j4 by —U'{x) — A. 

In order to prove Eq. dJ, let us consider separately the three regimes 
A £ [0,AI], A £ {Al,A^], and A £ {A^,oo). For A G [0,^5;] the result in 
Eq. (|14() is trivial, since Vad(^) = "^-i^) = '^+(^) = 0- the intermediate 
regime A G (^^,^2]) l^t us assume that the particle is initially located at 
a minimum of the potential U{x) + Ax (as we have mentioned before, the 
drift velocity does not depend on this particular initial condition). If we 
choose the time period of F(t) as T = 2Nt-^{A), with = 1, 2, 3, . . ., 
then, during the first half-period, Nt+{A), the particle stays trapped at 
the initial location. After the second half-period, 2Nt-^-{A), the particle 
arrives to a new minimum of U{x) + Ax separated from the initial location 
by a distance A^. During the third half-period, 3A^T+(yl), the particle stays 
trapped at that minimum, and so on. Thus, for a given A £ (^^,^2], all 
the frequencies 
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lead to the same value of the drift velocity 

Consequently, taking into account that liniTv^oo f^Tvl^) = 0, it follows that 

v^A) = hm V [A,u:n{A)] = (21) 

This proves Eq. H14() for A £ (A*, Ag], since in this regime V-{A) = 0. 

Finally, for A G {A2,oo), it is easy to prove that t^{A)/t^{A} is a 
continuous strictly increasing function of A which takes values in the in- 
terval (0,1). Let us assume first that for a given A E (^2,00) the ratio 
t^{A)/t-{A) is a rational number p/q, with p < q. Then, if we choose 
the time period of F{t) as T = 2NqT+{A) = 2NpT-{A), with iV = 1, 2, 
3, . . ., the particle travels a distance equal to N{q — p) every time period. 
Consequently, all the frequencies 

= m^) = N^) ^''^ 

lead to the same value of the drift velocity 

V [A,u;^{A)] = ^^^^^^^^ = - = ^ • (23) 

Taking into account once again that lim^r^oo ^n{A) = 0, it follows that 

vUA) = lim V [A, u;m{A)] = ^+(-^) + ^-(-^) , (24) 

Since any irrational number can be approximated as closely as desired by 
a sequence of rational numbers, Eq. (|14|) is also valid if tj^[A) /t-{A) is an 
irrational number. 

In Fig. IS] we represent the adiabatic drift velocity in Eq. (|14|) as a func- 
tion of the driving amplitude A. We have also plotted the numerical drift 
velocity obtained for u; = 0.01. Both results are indistinguishable within 
the resolution of the plot. 

Finally, from the above argument, another conclusion follows. As men- 
tioned before, the adiabatic value Uad(^) not only is reached in the limit 
a; ^ 0, but also for the finite frequencies given by Eqs. (fT^ and ((^ . This 
result becomes manifest in Figs. |21 and 01 as an infinite number of inter- 
sections between the horizontal line (adiabatic limit) and the drift velocity 
curve. These frequencies are at the core (see the red regions in Fig. [^J of 
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the above discussed fingers. Furthermore, they correspond to exact dis- 
placements of n spatial periods every driving period T. This fact, together 
with the interpretation provided by our analysis of the zero-current bands, 
lead us to infer that each finger corresponds to motion in which the particle 
advances a distance of about n spatial periods at those time intervals when 
F{t) = A. We have checked this fact in the numerical simulations. 

4. Conclusions 

To summarize, we have studied a one-dimensional deterministic ratchet 
model in the overdamped regime driven by a rectangular force. We have 
performed extensive numerical calculations, which allowed us to present 
a colored phase diagram showing the drift velocity as a function of a wide 
range of the driving parameter values A and lo. In addition, we have studied 
the conditions leading to bounded particle motion, obeying asymptotically 
periodic, bounded motion, x{t + T) = x{t). We have shown that the result- 
ing regions cover most of the phase space where there is no directed current. 
Moreover, we have provided an explicit analytical expression for the drift 
velocity in the adiabatic limit. Finally, the theoretical considerations are 
of great help to rationalize the finger-shaped structure shown in the diagram. 

We (DC, JC-P, AA, MM) acknowledge the support of the Ministerio de 
Educacion y Ciencia of Spain (FIS2005-02884) and the Junta de Andalucfa. 
DC also acknowledges the Ministerio de Educacion y Ciencia of Spain for 
a contract under the Juan de la Cierva program. The authors wish Profes- 
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personal happiness. 
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